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We  present  an  experiment  that  measured  the  behavior  of  a  shock  precursor  in  a  well-characterized,  thin  thermal  layer  above 
a  3000 “K  hot  surface.  A  laser  pulse  focused  onto  a  thin  solid  target  produced  the  shock,  eliminating  interference  from  explosion 
products.  Precursor  and  thermal  layer  characteristics  were  measured  using  spectroscopy,  interferometry,  and  dark-field 
shadowgraphy. 
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PRECURSOR-SHOCK  IN  A  WELL-CHARACTERIZED  THERMAL  LAYER 

A  thermal  shock  precursor  foms  when  a  shock  wave  propagating  through 
ambient  air  encounters  a  layer  of  hot  gas  whose  sound  speed  is  higher  than  the  sound 
speed  in  the  ambient.  Precursor  formation  begins  when  a  portion  of  the  incident  shock 
wave  refracts  into  the  hot  layer  creating  a  toe  wave  that  accelerates  ahead  of  the  main 
shock  and  increases  the  static  pressure  in  the  shocked  hot  gas  (Figure  1).  In  turn,  this 
increased  pressure  drives  a  shock  -  the  precursor  -  back  into  the  ambient  air.  When  the 
precursor  encounters  the  original  shock  wave,  a  reflected  shock  propagating  into  the  hot 
layer  is  formed.  The  interaction  of  this  reflected  shock  with  the  hot  gas  layer  generates 
vortical  flow  that  sweeps  additional  air  into  the  layer  and  reinforces  the  toe  wave  further 
to  complete  the  cycle.  Precursor  shocks,  which  were  first  observed  in  nuclear  explosions 
as  early  as  Trinity,'  are  thought  to  play  a  role  during  explosive  breakup  of  meteorites  in 
the  atmosphere,  high-velocity  impact  of  cosmic  bodies  on  planetary  surfaces,  and  in  the 
interaction  of  interplanetary  shocks  with  the  geomagnetic  tail.  Some  have  even 
suggested  harnessing  artificially  induced  precursors  to  alter  the  trajectories  of  small  but 
dangerous  asteroids  entering  the  earth’s  atmosphere.'* 


Very  few  experiments  have  measured  precursor  formation  in  well-characterized 
thermal  layers.’  Most  published  results  considered  precursors  near  the  interface  of  two 
gas  layers  stacked  above  a  plane  surface.  ’  The  gas  closer  to  the  surface  had  a  higher 
sound  speed  and,  therefore,  served  as  a  model  of  a  thermal  layer.  However,  the  two  gases 
in  some  experiments  were  kept  apart  by  a  plastic  sheet,  which  may  have  affected 
precursor  dynamics.  Also,  in  these  experiments  the  density  profile  of  the  higher-soimd- 
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speed  gas  either  decreased  away  from  the  surface  or  was  uniform,  a  situation  unlike  the 
density  profile  of  a  thermal  layer  in  which  the  surface  is  the  source  of  heat.  In  thermal 
layers  the  gas  density  is  lowest  near  the  surface  where  the  temperature  is  highest. 

In  this  paper,  we  describe  an  experiment  that  measmed  the  characteristics  of  a 
precursor  in  a  well-characterized,  thin  thermal  layer  produced  by  the  rapid  heating  of  a 
metal  foil.  The  brightness  temperature  of  the  hot  surface  was  measured  spectroscopically. 
Density  and  temperature  profiles  of  the  thermal  layer,  as  well  as  the  density  within  the 
shock  and  the  precursor,  were  determined  using  Abel-inverted  interferographs.  Dark-field 
shadowgraphs,  acquired  simultaneously  with  the  interferographs,  measured  the 
propagation  characteristics  of  the  shock  and  precursor. 

A  laser  pulse  focused  onto  a  thin  solid  target  produced  the  main  shock  wave  in 
our  experiment.  Shocks  produced  by  laser  pulses  were  proven  to  behave  like  ideal 
Taylor-Brode-Sedov  shocks  generated  by  more  conventional  means  such  as  explosions  or 

shock  tubes.’’*’  ®  However,  laser-generation  of  shocks  has  a  number  of  advantages. 
Pressures  achieved  with  laser-induced  shocks  can  be  much  higher  than  pressures 
achievable  with  either  chemical  explosives  or  shock  tubes.^  Moreover,  since  a  laser- 
explosion  has  little  mass,  the  shock  can  be  observed  for  a  long  time  without  being 
overrun  by  explosion  products. 
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The  experiment  took  place  inside  an  aluminum  box  of  ~  1 5-cm  length  per  side 
(Figrire  2).  A  hollow  cone  positioned  with  its  apex  protruding  into  the  box  was  attached 
to  one  of  the  box  walls.  The  cone’s  tip  was  cut  off  leaving  a  l-mm  diameter  opening  that 
was  closed  by  a  2.5-pm  thick  layer  of  mylar.  Inside  the  box  was  a  5  cm  x  5  cm  x  12  pm 
tantalum  sheet  oriented  with  its  surface  perpendicular  to  the  cone  axis  and  attached  in 
series  to  a  switched  780-pf  external  capacitor  charged  by  a  1.3-kV  power  supply.  The 
distance  between  the  cone  tip  and  the  tantalum  sheet  was  adjustable  from  zero  to  3  cm. 
Glass  windows  on  the  box’s  two  opposite  walls  allowed  viewing  and  access  to  optical 
diagnostics. 

To  produce  a  precursor  shock  the  box  was  first  placed  inside  a  larger  vacuum 
chamber,  both  the  box  and  the  chamber  were  evacuated,  and  the  box  was  back-filled  with 
one  atmosphere  of  argon  gas.  The  capacitor  was  then  discharged  through  the  tantalum 
foil,  heating  the  foil  and  a  thin  layer  of  argon  gas  above  it.  (The  parameters  of  the 
discharge  circuit  were  chosen  to  generate  a  fast-rising  current  that  heated  the  foil  quickly 
so  that  the  layer  of  gas  above  the  tantalum  remained  thin  and  laminar  for  the  duration  of 
the  experiment.  See  Figure  3a.)  One  millisecond  later  the  laser  was  fired  to  produce  a 
shock.  The  laser,  a  130-joule,  1.054-pm,  5-ns  duration  pulse,  was  focused  through  the 
vacuum  of  the  larger  chamber  onto  the  tip  of  the  cone  where  it  irradiated  and  rapidly 
heated  the  mylar  foil.  The  heated  mylar  expanded  rapidly,  much  like  the  products  of  a 
chemical  explosion,  forming  a  shock  in  the  argon  gas  that  propagated  into  the  box  and 
through  the  thermal  layer  above  the  heated  tantalum  sheet. 


3 


The  heated  tantalum  reached  brightness  temperatures  as  high  as  3000°K  (Figure 
3b).  The  temperature  and  its  time-evolution  were  determined  by  a  grey-body  fit  to  light 
emission  from  the  tantalum  surface,  spectrally  resolved  from  500-700  nm  and  recorded 
with  a  streak  camera.  The  light  emission  rose  to  a  maximum  in  40  psec,  together  with  the 
current,  and  remained  constant  within  13%  for  1  msec.  We  estimated  that  a  13%  emission 
variation  corresponds  to  a  temperature  variation  of  less  than  2%.  The  imaging  system 
was  calibrated  between  1500  and  2800°K  using  a  tungsten  lamp  with  emission 
characteristics  traceable  to  the  National  Bureau  of  Standards.  We  also  used  the  tungsten 
lamp  to  test  the  sensitivity  of  the  fitting  and  data  reduction  procedure,  from  which  we 
estimated  that,  despite  the  small  wavelength  range  measured,  temperatures  that  differ  by 
6%  could  be  distinguished.  Light  emission  from  the  tantalum  was  also  monitored  on 
each  calibration  shot  with  an  MRD  510  photodiode  sensitive  from  400  to  1 100  nm. 
Emission  in  this  visible  to  near-IR  band  was  correlated  with  foil  temperatures  determined 
spectroscopically,  and  the  photodiode  was  used  as  a  temperature  monitor  on  each  data 
shot. 


The  thermal  layer  and  the  shock  were  probed  with  a  0.53-pm,  0.35-ns  duration, 
and  5-cm  diameter,  50-mJ  laser-probe  that  illuminated  the  experiment  along  the  surface 
of  the  tantalum  sheet.  Laser  probe  light  passing  through  the  experiment  was  collected 
with  a  telescope,  imaged  using  dark-field-shadowgraphy  and  recorded  with  a  folded-wave 
interferometer. The  density  and  temperature  profiles  of  the  thermal  layer  are 
determined  from  the  probe’s  phase  shift  J given  by 
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-:;—  =  -\(Hz)-no)ds  =  (no-l)- - 1 

2^  Po  J 

where  «(,=  1.000282  is  the  index  of  refraction  of  argon  at  1  atm.  and  0°C,  L  is  the  length 
of  the  tantalmn  foil,  A  is  the  probe  wavelength  and  p(z)  and  Po  are  gas  densities  inside  and 
outside  the  hot  layer  respectively.  ’  ’  The  letter  z  designates  a  direction  perpendicular  to 

the  tantalum  foil  surface.  In  writing  this  expression  we  assumed  that  the  index  of 
refraction  n  is  proportional  to  density  with  negligible  explicit  dependence  on  temperature 
(i.e.  n  =  1  +  («o-l)p/po).  In  addition,  since  one  millisecond,  the  time  at  which  the  shock 
wave  was  initiated,  is  much  larger  than  the  microsecond  sound-transit  time  through  the 
hot  layer  (~  1  mm/319  m/sec),  the  pressure  in  the  hot  layer  should  be  equilibrated  with 
the  pressure  of  the  rest  of  the  gas.  Under  these  circumstances  the  temperature  T(z)  in  the 
layer  is  related  to  probe  phase  through 


(2) 
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where  Tq  is  the  temperature  of  the  unheated  gas. 

Figure  4  presents  the  temperature  profile  of  the  hot  layer  at  one  millisecond.  During  the 
30  psec  that  the  shock  wave  sweeps  over  the  tantalum  foil  (as  will  be  shown  below)  the 
layer  may  be  taken  to  be  stationary. 


Figure  5  through  Figure  7  show  development  of  the  shocks  in  our  experiment. 

The  principal  shock,  initiated  11.5  mm  above  the  surface,  was  first  observed  at  6.8  psec 
when  its  radius  was  17.8  mm:  In  the  subsequent  23.6  psec  it  propagated  to  32.6  mm.  The 
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precursor,  conversely,  was  2.5  mm  long  at  6.8  psec  and  developed  to  a  length  of  6.6  mm 
by  17.4  psec.  As  in  our  previous  work,’’*  ®  the  radius  r  of  the  principal  shock  scales  very 
well  as  the  2/5“*  power  of  time  t,  in  conformance  with  the  Taylor-Brode-Sedov  similarity 
solution  for  a  massless  point-source  explosion 

(3) 

V  Po  ) 

where  E  is  the  explosion  energy,  is  the  mass  density  of  the  ambient  gas,  and  a(y)  is  an 

adiabatic  index  dependent  correction  factor  related  to  the  efficiency  of  converting 
explosion  energy  to  shock  energy.'^  If  E  is  identified  with  the  laser  energy  then  an  a  of  2 
is  required  to  fit  the  absolute  r  vs.  t  values  to  the  data,  indicating  that  the  post-shock  y  is 
about  1 .2  in  the  experiment.'''  From  expressions'^  for  shock  pressure  P  and  shock 
velocity  r 


(4) 


P  = 


y-1 

2y 


(5) 


r 


5t 


where  c,  the  soimd  speed  in  argon  at  room  temperature,  is  330  m/sec,  we  infer  that  at  the 
first  observation  the  pressure  of  the  principal  shock  was  13-17  atmospheres,  for  a  pre¬ 
shock  y  between  5/3  and  1  respectively,  and  at  the  last  observation  2-3  atmospheres. 
More  precise  Hugoniot  calculations  using  the  SESAME  5172  equation-of-state  table  for 
argon  yield  pressures  very  close  to  those  obtained  from  equation  4  with  a  pre-shock  y  of 
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5/3:  At  6.84  |asec  they  give  a  pressure  of  13.3  atm,  at  21.3  psec  a  pressure  of  8.4  atm, 
and  at  30.4  psec  a  pressure  of  2.0  atm. 

During  the  observation  period  the  precmsor  height  defined  at  the  triple  point, 
(where  the  preeursor  shock,  the  main  shock,  and  the  reflected  shock  meet),  increased  at 
190  m/sec,  changing  from  2.1-mm  at  6.8  psec  to  6.6  mm  by  30.4  psec.  Interestingly,  the 
ratio  of  precursor  height  to  main  shock  radius  was  not  constant,  as  would  be  the  case  if 
the  precursor  propagation  were  self-similar.  Rather,  the  ratio  increased  in  an 
approximately  linear  fashion  from  0.1  at  6.8  psec  to  0.2  by  the  end  of  the  experiment. 


Theories  that  describe  the  formation  of  a  precursor  at  the  interface  of  two 
dissimilar  gases  of  density  Po  and  predict  that  the  angle  of  the  precursor  with  respect  to 


the  horizontal  (P  in  Figure  7)  is  given  as  sin(y^)  = 


1— It  is  interesting  to  apply  this 
Po 


expression  to  our  experiment  and  determine  where  the  analogous  densities  occur.  In  our 
experiment,  precursor  angle  is  24  degrees,  which  corresponds  to  a  density  ratio  of  0.17  in 
the  above  expression.  (The  angle  may  be  decreasing  slightly  from  26  to  23  degrees  in 
time.)  The  equivalent  temperature  ratio  in  our  experiment  would  be  6,  which  occurs  at 
about  1  mm  away  from  the  tantalum  surface  (Figure  4). 


The  density  profiles  inside  and  outside  the  shocks  are  determined  from 


interferometry  analyzed  using  the  inverse  Abel  transform  method.'®  The  technique  is  here 
applied  as  follows:  The  fringe  shift  associated  with  a  ray  traversing  heated  or  shocked  gas 
is  given  by  (Figure  8) 
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(6) 


S{r,z) 

In 


1 

-  ^(n{r,z)-n^{z))ds . 


Here  r  is  a  cylindrical  radius  extending  from  the  axis  of  symmetry.  If  z  is  above  the 
thermal  layer  then  (z)  =  is  index  of  refraction  of  unheated  argon.  But  if  z  is  within 
the  thermal  layer,  then  itf^z)  is  index  of  refraction  of  the  heated  but  imshocked  thermal 
layer  in  Figure  4.  (An  underlying  assumption  is  that  the  thermal  layer  is  reproducible 
from  shot  to  shot.)  Assuming  cylindrical  geometry,  and  noting  that  =r^  -x^ , 

sds  =  rdr ,  and  ds  =  rdr  /  -x^  (Figure  8),  the  integration  along  s  in  can  be  rewritten 

as 

d{r,z)  r  {n{r,z)-n^{z))2rdr 
2;r  * 

Using  the  Abel  transform  this  expression  can  be  inverted  to  yield 


(8) 


n(r,z)  =  no(2) - 

n 


dS(r,z)  dx 
2ndx  -r^ 
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which  can  be  numerically  solved  to  determine  the  index  of  refraction  n(r,z).  Furthermore, 
since 

(9)  w(r,  z)  =  1  +  («o  (z)  - 1)  , 

Po 

the  density  profile  can  be  determined  from 


pjr,  z)  _  ^  _  1  A  p  dS{r,  z)  dx 

Po  «o(2)-l^J,  27tdx 
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Determination  of  5(r,z)  requires  a  certain  amount  of  caution  and  judgment  since  the 
fringe  pattern  can  be  quite  complicated  as  shown  in  Fig.  5.  In  some  regions  of  the 
interferograph  it  is  difBcult  to  clearly  identify  and  associate  phases  with  fringe  lines. 

The  results  of  the  analysis  are  shown  in  Figs.  9  through  12.  Figures  9  and  10 
pertain  to  the  shots  with  a  thermal  layer  in  Figine  5,  while  Figs.  1 1  and  12  pertain  to  the 
shots  without  a  thermal  layer  in  Figure  6.  In  Fig.  9  we  show  grayscale  contoin  plots  of 
the  density  gradients  (left-hand  side)  and  the  density  p/po  (right-hand  side)  at  6.8, 10.2, 
17.4,  and  30.4  psec.  In  this  figure,  dark  gray  indicates  high  relative  densities,  while  light 
gray  indicates  low  relative  densities.  The  thermal  layer  is  analyzed  for  all  the  thermal 
shots  except  that  at  17.4  psec;  only  the  fringes  above  the  'wave-like'  discontinuity  in  the 
interferograph  at  17.4  psec  in  Fig.  5  could  be  clearly  identified  and  labeled. 

The  shadowgraphs  in  Figure  5  compare  well  to  the  reconstructed  density  gradient 
contours  in  Figure  9.  The  major  shock  structure  around  the  triple  point  is  obtained  from 
the  Abel  analysis,  as  well  as  some  of  the  minor  structures  (e.g.,  at  6.8  psec).  One 
problem,  somewhat  evident  in  Figure  9,  especially  at  6.8  and  30.4  psec,  is  a  noisy  hash 
seen  most  clearly  in  the  upper  right  portions  of  the  density  gradient  plots.  This  arises 
because  of  the  nature  of  the  data  and  the  inverse  Abel  transform.  The  integrand  of  the 
inverse  Abel  transform  is  proportional  to  d5ldx\  when  the  data  is  noisy  this  can  lead  to 
large  variations  in  the  slope  of  5  and  hence,  noisy  results.  Some  data  smoothing  is  used 
to  mitigate  this  effect. 
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The  density  contours  of  p/p^  in  Figure  9  show  that  the  gas  undergoes  a 
compression  at  the  shocks  and  then  a  subsequent  rarefaction  (this  latter  point  is  most 
obvious  at  6.8  psec).  The  amount  of  compression  appears  to  be  a  function  of  relative 
position  to  the  triple  point  where  the  density  enhancement  is  largest.  The  actual  relative 
density  profiles  p/p^  as  a  function  of  radius  and  height  above  the  plate  (z  direction)  are 
shown  in  Figure  10  at  6.8, 10.2, 17.4,  and  30.4  psec.  The  z-direction  in  this  figure  also 
measures  the  relative  density;  the  spacing  between  tick  marks  is  unity.  For  example, 
consider  the  topmost  density  profile  in  Figure  10  at  6.8  psec,  which  corresponds  to  the 
top  of  the  density  contour  in  Figure  9.  The  gas  is  undisturbed  near  the  left  boimdaiy  so 
that  p/po=  1.  At  the  shock,  the  relative  density  rapidly  increases  to  slightly  more  than  one 
unit,  so  that  p  /po  >  2.  Continuing  downstream  of  the  shock  (to  the  right  in  the  figure),  the 
density  p/pg  decreases  below  imity.  The  maximum  density  jump  obtained  from  the 
analysis  is  p/pg  s  2.9  and  occurs  at  6.8  psec.  Hugoniot  calculations,  mentioned  earlier, 
give  p/po  =  3.1  at  6.8  psec,  p/pg  =  2.7  at  21.4  psec,  and  p/pg  =  1.5  at  30.4  psec,  which  is 
the  same  result  one  would  get  fi'om  the  expression'^ 


(11) 


P_ 

Po 


y  +  \ 

y-\ 


[1  + 


2 

y-\ 


withy  =  5/3. 

Figures  1 1  and  12  are  similar  to  Figs.  9  and  10,  but  are  for  shots  without  a  thermal 
layer  (see  Fig.  6).  Because  there  is  no  thermal  layer,  the  shock  structures  are  not  very 
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complex.  Again,  the  density  gradient  contours  in  Fig.  11  are  to  be  compared  with  the 
shadowgraph  in  Fig.  6.  The  agreement  between  the  two  is  good,  especially  at  17.4  psec. 
At  30.4  psec  the  shadowgraph  (Fig.  6)  shows  an  almost  horizontal  discontinuity 
downstream  of  the  triple  point.  Although  there  is  some  indication  of  this  in  Fig.  1 1,  it  is 
very  weak  and  comparable  to  noise  levels.  The  density  contours  in  Fig.  1 1  also  show  a 
strong  compression  at  the  shock  followed  by  a  rarefaction  of  the  gas.  At  17.4  psec  there 
is  a  compression  and  rarefaction  associated  with  the  first  shock  stracture,  followed  by  a 
weaker  compression  and  rarefaction  with  the  second  shock  stmcture.  These  features  are 
also  shown  in  Fig.  12  which  is  a  plot  of  the  relative  density  profile  p/po  as  a  fimction  of 
radius  and  height.  The  maximum  density  compression  is  p/po  s  2.7  at  17.4  psec. 

The  authors  would  like  to  thank  Mr.  Levi  Daniels  for  his  technical  support  and 
Ms.  Naomi  Frankel  for  analyzing  some  of  the  data.  We  are  also  grateful  to  Dr.  Gary 
Schneyer  of  Maxwell  Inc.  and  Mr.  Emil  Braun  of  NIST  for  helpful  discussions  and 
suggestions.  This  work  was  supported  by  the  Defense  Special  Weapons  Agency. 
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Figure  1 


toe 


Schematic  of  precursor  shock  formation.  Shock  propagating  through  air  refracts  in  a 
thermal  layer  creating  a  toe  wave  that  accelerates  ahead  of  the  main  shock.  Increased 
pressure  behind  the  toe  drives  a  precursor,  which  upon  reflection  from  the  main  shock 
creates  a  reflected  wave  that  propagates  back  towards  the  thermal  layer.  The  interaction 
of  this  wave  with  the  thermal  layer  generates  vortical  flow  that  sweeps  additional  air  into 
the  thermal  layer  and  further  reinforces  the  toe. 
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Figure  2 


Tantalum  foil 
5cm  X  5cm  X  12  ^m 


to  scope 


laser 


Setup  of  experiment 


Time  (msec) 


(a)  Profiles  of  the  tantalum  foil  current  and  the  light  emission  from  its  surface,  (b) 
Temperature  of  the  tantalum  foil  vs.  charging  voltage. 


Figure  4 


Temperature  profile  of  Argon  gas  as  a  function  of  distance  from  the  tantalum  foil  at  one 
millisecond.  The  profile  was  determined  using  interferometry  (inset).  When  the 
temperature  is  large  (and  density  small)  small  phase  errors  translate  to  large  temperature 
errors. 
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Dark-field  shadowgraphs  and  interferographs  showing  time-evolution  of  a  precursor 
shock.  Ground  range  on  the  shadowgraphs  is  measured  in  centimeters.  The  laser  energy 
was  138,  104,  87,  and  137  joules  respectively.  The  blast  was  initiated  1 1.5  mm  above  the 


Figure  6 


17  lasec 

no  thermal  layer 


30  i^sec 

no  thermal  layer 


Shocks  at  17  ^sec  and  30  ^sec  over  a  surface  with  no  thermal  layer.  Ground  range  on  the 
shadowgraphs  is  measiued  in  centimeters.  The  laser  energy  was  83  joules,  and  135  joules 
respectively.  The  blast  was  initiated  11.5  mm  above  the  surface. 
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Figure  9 


Grayscale  contour  plots  of  the  density  gradients  (left  hand  side)  and  the  density  ratio  p/po 
(right  hand  side)  at  6.8, 10.2, 17.4,  and  30.4  psec  for  the  shots  shown  in  Figure  5.  The 
dark  gray  indicates  high  relative  densities,  while  the  light  gray  indicates  low  relative 
densities.  The  density  gradient  plots  are  to  be  compared  to  the  shadowgraphs  near  the 
triple  point  in  Figure  5.  20 


Figure  11 


Grayscale  contour  plots  of  the  density  gradients  (left  hand  side)  and  the  density  ratio  p/po 
(right  hand  side)  at  17  and  30  psec  for  shots  without  a  thermal  layer  shown  in  Figure  6. 
The  dark  gray  indicates  steep  gradients  and  high  relative  densities,  while  the  light  gray 
indicates  weak  density  gradients  and  low  relative  densities.  The  density  gradient  plots  are 
to  be  compared  to  the  shadowgraphs  near  the  triple  point  in  Figure  6. 
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